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Abstract. We investigate the existence of several (anti-)ferromagnetic phases in the 
diluted ferromagnetic Kondo-lattice model, i.e. ferromagnetic coupling of local moment 
and electron spin. To do this we use a coherent potential approximation (CPA) with a 
dynamical alloy analogy. For the CPA we need effective potentials, which we get first 
from a mean-field approximation. To improve this treatment we use in the next step a 
more appropriate moment conserving decoupling approach and compare both methods. 
The different magnetic phases are modelled by defining two magnetic sublattices. As a 
result we present zero-temperature phase diagrams according to the important model 
parameters and different dilutions. 
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1. Introduction 

In solid state physics magnetism is one of the most discussed phenomena. Especially 
for spintronic devices, which use in addition to the charge the spin of the electrons, it 
is necessary to understand the magnetism of highly correlated materials. One class of 
materials that could be used for spintronics are the diluted magnetic semiconductors 
(DMS) ( [H [21 El H] ) . These consist of semiconductors as a host material in which small 
amounts of magnetically active atoms are deposed. The ambition is that the resulting 
material has the properties of a semiconductor and is magnetic, too. But the mixture 
of materials can also destroy one or both qualities. 

To describe DMS in theory several approaches are used. Much work is done based on ab 
initio calculations El [TJ El E] - Those calculations can give material specific information 
and therefore show the chemical origin of differences between various semiconductors. 
In addition the exchange mechanism is usually described by a model Hamiltonian. 
It lies in the nature of a model to simplify the actually appearing physical processes 
and to concentrate on the most important mechanism occurring in a certain material. 
Thus diverse models are used to describe different DMS. These include, e.g., the double 
exchange [§], superexchange [TO] and/or the Kondo-lattice model (KLM) [TT1 [T2l [T3] . 
Since the KLM is widely used for the description of DMS, we want to discuss its possible 
magnetic ground states. Doing this we want to compare the internal energy of different 
states. This, of course, implies a predetermination of those states. Besides the fully 
ordered ferromagnetic and the completely disordered paramagnetic state, regarding the 
local moments, there can be, e.g., several antiferromagnetic states, spin-canted ones, 
spirals, spin-density waves or even spin-glas phases. All of them were found to be 
ground states in specific models or parameter regions [HI [151 IB]- Since a complete 
description of those states is not possible in our model we will focus on five important 
types of possible ground states. This should yet show some general trends of the phase 
diagram of the diluted KLM. Zero temperature magnetic phase diagrams were already 
calculated for the concentrated KLM ([TTl HH HU [20]), but, to our knowledge, not for 
the diluted one. 

This work is organized as follows: First we show which magnetic phases are investigated 
in this paper and how they are described by a definition of magnetic sub-lattices. Then 
we use a equation of motion approach to get approximate solutions of the many-body 
problem of the Kondo-lattice model. This gives us effective potentials which we can 
include in a coherent potential approximation. After that we are able to calculate the 
internal energy of the single magnetic phases and present phase diagrams for different 
concentrations. The origin of those phases are discussed by means of quasi-particle 
density of states. Finally we summarize this work and give a short outlook to future 
work. 
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2. Model and Theory 

We want to investigate which type of magnetic order is present in diluted magnetic 
systems at zero temperature. In our work the magnetic component is described by a 
KLM with a positive coupling J, commonly known as ferromagnetic Kondo-lattice model 
(FKLM) ( [2T1 122] ) . although the notion may be considered as somewht misleading. Since 
the carriers in this model study are electrons, this positive sign of J is equivalent to 
a negative one in the hole picture |23J which is usually assumed for DMS. The FKLM 
is characterized by an interaction between itinerant electrons of spin cr^ and localized 
moments S». Its Hamiltonian reads 

H = H s + H sd = t^2 c ^ c 'J a ~ J ' <Ji W 

(ij)o i 

Due to the interaction part H S( i it is energetically favorable when the electron's and the 
local moment's spin on the same site R, align parallel for positive J. The next- neighbor 
hopping t of the itinerant electrons leads to an indirect coupling of the local moments 
on different sites. For our purposes it is better to represent the spin of the electrons by 
their creators/annihilators via <Tj = | c ta^^' c io'^ where cr is the vector of the 
three Pauli matrices. When we write down the scalar product Sj • <Xj explicitly we get 

H — t Cio C jo- ~ 5^ X i { Z o-Si cf a Ci a + SI C^i^Cia) . (2) 

(ij)cr i<r 

S? is the ^-component of the Sj and Sf = Sf + z a iSf , with z a = 5 a ^ — 8 C ±, are the ladder 
operators. To take into account disorder we use a site dependent parameter Xi = 0, 1 
which turns the interaction on or off. The sum over these parameters is fixed by the 
given concentration c = 4 £V of the magnetic atoms. 

To include antiferromagnetic phases we have to divide the lattice into sub-lattices. 
In this work we only consider (anti-)ferromagnetic phases which can be divided into 
two magnetic sublattices (cf. figure CD). Within each sublattice we can define a 
sub-lattice magnetization, where the magnetizations of the two sub-lattices for the 
antiferromagnetic phases are antiparallel to each other. We assume a ferromagnetically 
saturated sub- lattice magnetization at T = 0. Even though it is known that these 
Neel states are not the ground states of the model, due to spin fluctuations, they are 
energetically close to the real ground states, especially for large spins. To take the 
sub-lattice formulation into account we change the Hamiltonian (J5J) to 

H = ^ ] ^cf aa Cjp a — — ^ ] Xi a {Zo-S^C^^iaa + S ia C i _ aa Ci aQ ) . (3) 
(i,j)craf3 icra 

The lattice sites Rj, Rj are now within a magnetic sublattice and are shifted by a 
primitive translation vector r a according to the specific sub-lattice. Since we yet have 
next-neighbor hopping in the chemical lattice the hopping within each sub-lattice is 
constricted. For example the next-neighbors of the magnetic and the chemical lattice 
in the A- type phase are only the same in, e.g., xy-planes. Thus the intra-lattice hopping 
if? is only two-dimensional. The hopping in ^-direction leads to a hybridization of the 
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Figure 1. (Color online) Magnetic phases considered in this work. 



two sub-lattices. Fourier-transformation of the hoppings yields the dispersion of the free 
system and can be written for the A-type AFM as 
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(cos(ak x ) + cos(aky)) 



W 

K (k) = — cos(a/c 2 ) 



(4) 
(5) 



where W = 12t is the free bandwidth of the bulk system and a the lattice constant. For 
the C-type the intra-lattice hopping is one- dimensional and in the G-type there is no 
hopping within the magnetic lattice. 

The ferromagnetic phase needs no definition of sub-lattices, of course. But this 
could be also done artificially by choosing one of the lattice decompositions of the 
antiferromagnetic phases. The sub-lattice magnetizations are now parallel ((Sj Q ) = 
(Si a >)) to each other. No matter which decomposition we choose, we get the same 
results as for the non-decomposed case. This is a strong confirmation of the validity of 
the sub-lattice method. 

2.1. Internal energy 

The internal energy U of the FKLM at T = is given by the ground-state expectation 
value of the Hamiltonian ([3]). A straightforward calculation using the spectral theorem 
then shows, that U can be calculated from the retarded local electronic-GF (G^f a (E) = 



U=(H) = J2 T Ef-(E) (--ImGTAE)) dE, 



(6) 



where /_ (E) denotes the Fermi function. Until now this is an exact formula but since ([3]) 
defines a complicated many-body problem approximations in calculating the electronic- 
GF have to be accepted. Therefore we come now to the discussion of this point. 



2.2. Treatment of the many-body problem 

Our starting point is the equation of motion (EOM) of the electronic GF: 

J 



E W - tT")Gt a \E) = 5f - x t - (zJ%:\E) + F^f(E)) , (7) 

la" 



with Ising-GF 



(8) 
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and spin-flip-GF (SF-GF) 

K a k -:\m = ((srjc k ^ a ,c+ aa ,)) . (9) 

The main idea for construction of the phase- diagram is to assume certain configurations 
for the local moment system and then to compare their respective internal energies. 
Under the assumption of saturation the Ising-GF can be decoupled directly: 

c taa')) -> z <* SG iji{ E )- In PS] we nave shown, that the influence of the 
SF-GF on the ground state phase-diagram is small in the case of a system without 
dilution, at least if we only consider phases which are ferromagnetically saturated in 
each sub-lattice. In a first attempt to solve (JTj) we thus neglect the SF-GF arriving at: 

E ((£ + x lZa z^S)6T" - C") Gff MF) (E) = Sf. (10) 
We will call this solution the T = mean-field (MF) solution. 

To investigate the influences of spin-flip processes we will also perform a moment 
conserving decoupling approach (MCDA) as in Ref . [21] . One advantage of this method 
is, that we are not restricted to ferromagnetic saturation within the sub-lattices. Thus 
we can also specify the paramagnetic phase by setting (Sf) = 0. We formally define a 
self-energy in (JTJ) which reads as 

E W - tr")Gt a \E) = S$ + (([c_, H sd U c^J) 

la" 

= C' + * E M Sf(E)Gt a '(E). (11) 

la" 

The self-energy will play the role of a random potential later on, but we will first focus 
on a concentrated lattice (xi = IVi) to derive the formula for this self-energy. We assume 
M£"" (E) to be local which will not give rise to any contradiction in our treatment [24J. 
Thus we set M^°" (E) = M aa (E)5°f l " and have the simplified EOM of the concentrated 
lattice 

E (( E - M ^E))5T" ~ tT") G?;/(E) = 5?f . (12) 

la" 

When we compare ((7|) and ( TTTT) we find a relation of the self-energy and the higher 
Green's functions 

M acr (E)G^(E) = - J - (zJ^'(E) + F^\E)) 

= ~ { (zASL)G??M + z a T%:'(E) + FggiE)) (13) 

with the reduced Ising function Tf^' (E) = '(E) - (S z ia )G^(E) = ((5S*c i(ra] c+ aa ,)). 
Furthermore we read off from (TTTT) a formal correlation between self-energy and 
commutator 

[c iaa ,H sd ] — y M aa (E)c iaa , (14) 
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being rigorous, of course, only within the respective Green's functions. Now we write 
down the EOM for the higher GFs 

E( M ^''-c'')^ w (^) = 

la" 

(([S£, H sd U ka ^- c+ + ((Sr°[ Cka _ a , H sd U c + J) , (15) 

E(^ a "-c") nt a \E) = 

la" 

(([8St a , H sd }.c kaa , c+,„» + (m a [c katT , H sd }_- c + (16) 

We will perform different treatments of the non-local (i ^ k) and the local (i = k) terms 
in these equations. Let us start with the non-local ones. It has been shown in [23] that 
the commutators containing a spin operator are connected to magnon energies. Those 
are typically several orders of magnitude smaller than the electron energy. Thus we 
neglect these commutators. For the other ones we assume that we can replace them by 
the formal expression ( I14p . 

This procedure is however inappropriate for the local terms due to the intra-atomic 
interaction between electron and localized spin on the same lattice site. In this case we 
consider the commutators explicitly. Equations (|T5|) and ( IT6|) now read 

£ ((£ - M„ aa {E)Kf ~ t a k f)F™: a \E) = 



la" 



E 

la" 



8g (-M_ aa {E)F™:\E) + (([SrJ Ci ^ a , H sd }_; c+ , (17) 



8g (-MnWTtg? '(E) + ({[6S* a c iaa ,H sd Uc+ aa/ ) 

where we avoided double counting of the self-energy for the local terms by subtracting 
it on the right-hand-side. Multiplying ^ ka G^ k \ a (E) to the according equation leads 
with (CCD to 

F%?{E) = G™ a {E) (-M_ aa (E)F^'(E) + (([^, HJ.; c+ Q ,») (19) 

T™«'{E) = GTAE) (-M aa {E)T™« '(E) + (([5S* a c i<7a , H sd }_; c+ Q ,») . (20) 

The remaining commutators create various Green's functions. Additionally to the known 
ones we get four higher GF 

Euijcr (E) — ({S ia S ia Ci- aa ; Cj ffa if), 
F«"T\E) = ((5(S-°S^)c iaa ;ct a ,)), 



paa'(3)/p\ // q— a . + 

^iiiija \^ ) ~ \Wia n icraCi-aa] C j a a' 



Eiiiija (E) — ((SiarLi-a-aCiva; Cj aa ,)). 

In ferromagnetic saturation we can express two of them exactly 



W(1) ' • / - ! 1 



(E)=[S-- + -z a )F^(E), (21) 
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F^\E) = 0. (22) 



ill] a 



This cannot be done for F°^^'^ (E) but we know that in the limiting case of an empty 
conduction band both functions vanish and for a completely filled band we find 

F°^(E) = F%? '(E) (23) 

Fg^(E) = T^'(E) + [St)G^{E) . (24) 

Thus we make the ansatz that these functions can be expressed for all electron densities 
n and magnetizations (S z ) by the linear combinations 

Fuu£\E) = o&G${E) + ^F^'(E), v = 1, 3 (25) 

Fg^\E) = a^G^(E)+/3i^:'(E), a* = 2,4. (26) 

We fix the coefficients by the requirement that the zeroth spectral moments of the 
functions are conserved 



M^(F^\E)) = a^M^{Gtf c {E))+^M^{X^'{E)), v = 1,2,3,4 (27) 



with X = F, r. To determine pad we need additional conditions. For the first two 
functions we get these from the next moments 

M W (F^(E)) = *ttMW(G^(E)) + f}gM*\X%? iE)), u = 1,2 . (28) 

Furthermore we can calculate expection values via the spectral theorem from the 
different GF, which should also fullfill ( 125]) and f )26|) . Replacing the GF by the according 
expectation values yields 

F '■ (£>i a c iaa n io!aCia-a) = = Ot a J (Tli aa ) + P a }{S ia C iaa Ci a - a ^ (29) 

= a { a 4 l(n ia . a ) + &l a (5S z ia n ia _ a ) . (30) 

The first expectation value vanishes due to the combination of the electronic operators 
and the other is spin invariant. Equation ( !3~0|) seems to lead to no solution, due to 
underdetermination of the parameters. But ( 1271) already yields oial = (Sf a ni a _ a }, since 
the zero-moment of the spin-flip function F^f (E) disappears. Thus, as the only general 
choice for the other parameter, remains (3a} = (nj a _ CT ). The explicit terms of the other 



coefficients can be found in Appendix A With these conditions we can trace back all 
appearing Green's functions to a set of only three GF. Now we get a linear system of 
equations from ( T19T201) as 

•M SE'ff? ) = -io^(E) I | . (3d 




Tf^(E) J 2 
with the coefficient matrix 

I _[ ^- + Gf i c L a (E)(M a - (T (E) + ^C Fa a) ^B Faa Gi£ a (E) 

aa \ jC raa GT a (E) l + GZ(E)(M aa (E) + jB raa ) 
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The coefficients (A Faa , B Faa , . . .) are explained in Appendix A Within the matrix 



there are only local terms which are not site dependent when we assume translational 
invariance within the sub-lattice. Therefore we can write the solutions of (EH) as 



F?™'(E) = F?(E)G?g(E) , (32) 
Tl a 3 :\E) = T«{E)G«£{E) . (33) 

The only dependence on the lattice indices is in G?£(E). Comparing this with (fI3|) 
yields the self-energy 

M aa (E) = - 3 - {z a (S z ia ) + z a r:(E) + F^E)) , (34) 

which is indeed local. In an ordered system we would put the self-energy into formula 
(Fl2|) and had a closed set of equations. For a disordered system we will change this 
procedure a bit as described in the next section. 

2. 3. Inclusion of disorder - CPA 

The formulas derived in the previous chapter still depend on the concrete placement of 
the local moments. For a macroscopic large system the physical properties should be 
rather described by averaged quantities as long as the system is far from critical behavior 
near phase transitions. Therefore our goal is to find the configurationally averaged GF: 
(G^(E)). The problem is solved if we are able to find the self-energy T^ (T (E) of an 
effective medium with the according symmetry of the considered type of magnetism. Its 
general formula is 

Gff a \E) = (Ggm = ((E - E£)«W - e^'(q))" 1 . (35) 

V / aa' 



Let us inspect the MF-GF ffTO]) first. The disorder problem is defined by a random 
change of the local potential r] a for magnetic (M) and non-magnetic (NM) sites: 
■q^ a = —z a z a ^S; = with concentration c and 1 — c respectively. This is analogous 
to a binary alloy (for each direction of electron spin) with local disorder. We have applied 
the coherent potential approximation (CPA)[25J for calculation of the effective medium 
self-energy. Since we are only interested in the diagonal parts of the GF G^ a ,aa (E) and 
because of the locality of CPA self-energies, we can use a standard derivation of the CPA. 
One obtains the following (self-consistent) equation of a local self-energy approximation 
for the effective medium: 

M _ ycS NM _ yeS 

= c i^z a ° +(l-c) — aa ff ■ (36) 

i - « - sfjGfr* i - (vt - ^i)Gfr 

Due to (1351) the effective self-energies Yf^ a {E) of the two sub- lattices are not independent 
of each other. 

To improve the CPA result we now choose the self-energy of the MCDA M aa (E) = 
rj a(T (E) in ( 1341) as an effective random potential. This potential is energy dependent, 
that is why we call this a dynamical alloy analogy (DAA, [26| 127]). But we have to keep 
in mind that self-energy was derived for a concentrated system. In a diluted system the 
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Figure 2. The self-consistency cycle of the combined MCDA/CPA treatment. The 
magnetic sub-system (index " M" ) and the whole diluted system (index " eff" ) are not 
calculated separately, but in a combined loop. 



magnetic subsystem hybridizes with the whole system. 

Let us assume that we have found a solution for the self-energy Yff(E) of the diluted 
system. Now the Green's functions of the subsystems can be derived via the projection 

EE] 



G-£ a (E) = na v ' = , v = M, NM . (37) 

^-{^{E)-^f a {E))Gfr{E) 



We choose as an approximation that the GF of the magnetic subsystem G^ a (E) 
can be treated as an effective concentrated system which has the same EOM as in a 
real concentrated system. Thus the formal structure remains the same and we just 
have to replace G^ ±a {E) in (ST]) by G™ ±a±a {E). The self-energy and the magnetic 
sub-system's higher Green's functions {T^°^ a ' ' (E), F^ aa ' (E)) are now functions of the 
projected Green's functions G^ a a a {E). To take into account the hybridization effects 
with the complete system we do not put the self-energy M aa (E) directly into formula 
(TTTj) to get G}j£ a (E) as we would do it in a concentrated lattice. Instead, we use it 
to get the effective self-energy Y^ a (E) defined in Eq. (1361) . With S^.(£ l ) we calculate 
G^ ia (E) and the projected Green's function ( 1371) and the whole set of equations is 
closed. This self-consistency cycle is depicted in figure [2j 



3. Results and Discussion 

3.1. Mean- field results without dilution 

The internal energy of the FKLM at T = is given as an integral ([6]) over the product of 
(sub-lattice) quasi-particle density of states (QDOS) times energy up to Fermi-energy. 
For understanding the resulting phase-diagrams it is therefore useful to have a closer 
look at the QDOS first. In figure [3]the sub-lattice MF-QDOS is shown for the different 
magnetic phases investigated. The underlying full lattice is of simple cubic type with 
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a-AFM c-AFM g-AFM 



Figure 3. Sub-lattice quasi particle density of states (QDOS) of up and down electrons 
obtained from the MF-GF (fT0|) for two values of local coupling J shown for different 



magnetic configurations. 
eV. 



Parameters: S 



and free electron bandwidth: W = 1.0 



nearest neighbor hopping t chosen such that the bandwidth W is equal to W — 1 eV 
in the case of free electrons ( J = eV). The local magnetic moment is equal to S — |. 
We have plotted the up and down-electron spectrum separately for two different values 
of J = 0.1/1.0 eV. The exchange splitting A ex = JS eV of up and down-band is clearly 
visible. The decisive difference between the phases for nonzero values of J is bandwidth 
reduction from ferromagnetic over a, c to g-AFM phase. The reason for this behavior 
becomes clear by looking at the magnetic lattices shown in figure [U In the ferromagnetic 
case an (up-) electron can move freely in all 3 directions of space without paying any 
additional potential energy. In A-type anti-ferromagnetic phase the electron can still 
move freely within a plane but when moving in the direction perpendicular to the plane 
it needs to overcome an energy-barrier A ex . Hence the QDOS for large values of J 
resembles the form of 2D tight-binding dispersion. The bandwidth is reduced due to 
the confinement of the electrons. In the c-AFM phase the electron can only move freely 
along one direction and the QDOS becomes effectively one dimensional. Finally in the 
G-type phase the electron in the large J limit is quasi localized and the bandwidth gets 
very small. This bandwidth-effect is the main reason for the phase diagram obtained as 
shown in figure HI 

For larger J (J > 0) a typical sequence appears: for low band-fillings n the system is 
always ferromagnetic and, with increasing n, it becomes A-type then C-type and finally 
G-type anti-ferromagnetic. This behavior is understood easily by looking at the formula 
for the internal energy ([6]) and the MF-QDOS in figure [3j Because of the bandwidth- 
effect discussed already, the band-edge of the ferromagnetic state is always lowest in 
energy and will give therefore the lowest internal energy for small band-occupation. 
But since the upper edges of the low energy spin-bands decrease in the order fm, a-, c-, 
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0.2 0.4 n 0.6 0.8 1 



Figure 4. Phase diagram obtained from the MF-GF (fTU)) without dilution regarding 
coupling J and electron density n. Parameters: S — | and free electron bandwidth: 
W = 1.0 eV. 

g-AFM, the antiferromagnetic phases become eventually lower in energy for increasing 
band-filling. Note that for n = 1 the lower spin sub-band is completely filled. Therefore 
the bandwidth-effect is the main reason for the observed order of phases with increasing 
n. A very interesting feature can be found in the region: J = 0.2 . . . 0.3. In this 
region the ferromagnetic phase is directly followed by the c-AFM phase for increasing 
n although the a-AFM phase has a larger bandwidth than the c-AFM phase. This can 
be explained by the two-peak structure of the c-AFM-QDOS. Due to the first peak at 
low energies these energies are much more weighted than in the a-AFM case and the 
c-AFM phase will become lower in energy than the a-AFM phase. Since the reduction 
of bandwidth of the anti-ferromagnetic phases compared to the ferromagnetic phase 
is more pronounced for larger values of J the ferromagnetic region is growing in this 
direction. 

The paramagnetic phase does not appear for any finite J. Due to the down-shift of the 
up-spectrum of the ferromagnetic sub-lattice its internal energy will always be lower. 
This will change when we include correlation effects within our MCDA treatment of the 
FKLM. 

3.2. MCDA results without dilution 

All effects that have been discussed for the mean-field results are also present in the 
MCDA. It means that the differences of the bandwidths and the peak structures of 
the single phases remain the dominant origin of differences in the internal energy. But 
mainly due to the spinflip GF, which vanishes for the mean-field approximation, the 
QDOS changes. Thus the phase diagram should be changed, too. Figure [5] shows the 
QDOS in the ferromagnetic phase of a concentrated lattice. The main change compared 
to the MF-QDOS is the appearance of spin-down states at the energies of the spin- 
up band. These states can be connected to magnon emission processes. It turns out 
that the spectral weight of this part of the spin-down QDOS is heavily affected by 
the spin quantum number of the local moments. The smaller the spin the stronger 
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E[eV] 



Figure 5. QDOS of the MCDA-GF in the concentrated system for various spins S 
(JS = leV) in the ferromagnetic phase. The smaller the spin the larger is the amount 
of spin scattering states in the spin-down QDOS at E w — iS. The upper sub-band 
at E r; ^(5 + 1) gets to lower energy for increasing spin and constant JS. For S — > oc 
the scattering states vanish and the mean-field QDOS is reproduced. Parameters: 
W = leV , n — 0.05 




Figure 6. Phase diagram obtained from the MCDA-GF (fTOj) without dilution. 
Parameters: S = § and free electron bandwidth: W — 1.0 eV. 



are the influences of the scattering states. For very large spins the mean-field picture 
is reproduced (classical limit). But as it can be seen in figure E] the phase diagram 
changes for spin S — f, especially at stronger couplings J. Actually we get a big 
region of the paramagnetic phase, which was never present in the mean-field diagrams. 
This is plausible because in the paramagnetic regime the spin-flip processes should be 
important, which are not neglected by the MCDA. Thus the energy of this phase is 
much lower than for the mean-field approximation and can be of the same order as in 
the (A)FM phases. This broadening of the PM phase with increasing quantum character, 
was found in Ref. [29], too. 

At low couplings the MF and MCDA phase diagrams are similar again, also for S = |. 
This can be easily understood if we arrange the MCDA self-energy (|34j) in powers of J 
yielding: 



MAE) = ~z a (S?) + ^M' a (E) . 



(38) 
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Figure 7. MF-QDOS of the ferromagnetic (left) and A-type anti-ferromagnetic (right) 
configuration for various concentrations c of local moments. Parameters: J — 1.0 eV, 
S 



| and free electron bandwidth: W — 1.0 eV. 



We see that the self-energy splits into a MF and a "many-body" part M'^E). Since the 
MF term varies linearly with J it becomes dominant at small couplings. 

3.3. Results with dilution 

After we made clear the differences and similarities of the MF and MCDA in an undiluted 
lattice we want to discuss the effects of dilution. These are in principle the same for 
both methods. In figure [7] we show the MF-QDOS for the fm and a-AFM configuration 
at various concentrations c of the local moment system. With increasing dilution 
the spectral weight of the magnetic sub-bands decreases (proportional to c) and an 
uncorrelated band appears around the band-center of gravity T = eV. The shape 
of the magnetic sub-bands will still resemble the form of the undiluted case (for not 
too small c) whereas the appearing non-correlated band is more or less the same for all 
magnetic configurations and only depends on the dispersion of the underlying chemical 
lattice. For large exchange splitting J the magnetic sub-bands are well separated 
from the non- magnetic band and we can define their effective filling n eG — -. If the 
magnetic sub-system would be independent of the total system the phase diagrams for 
all concentrations c would be the same according to n eS . But due to hybridization of the 
magnetic sub-system with the total system the antiferromagnetic phases are remarkably 
suppressed for n cS < 1 (Figured]). This is best recognizable for the G-type AFM. For 
low concentrations (here c = 0.4) this phase is almost vanishing for ^ < 1 although it 
was present in the undiluted system for n > 0.9. 

For fillings larger than n = c the non-magnetic band starts to be filled and the phase 
will not change any more. Thus for n > c only the g-AFM phase exists due to a 
completely filled lower correlated band corresponding to an effective half-filling at the 
correlated spectrum. From figures H] and [6] we know g-AFM is just the existent phase 
for an undiluted lattice at n = 1. 

As seen in the case of an undiluted lattice the differences of the MF approximation and 
the MCDA are largest for strong couplings. That is also true in the diluted systems. 
But again, as for c = 1, this is more prominent for low spins. Actually there are almost 
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Figure 8. QDOS of the MCDA-GF Gf a aa (E) (bold black line), the (non-) magnetic 
sub-system's GF G^ a (E) (red line), G™' aa {E) (green dashed line) in a diluted 
system (ferromagnetic phase). For better comparing the sub-system's GF are 
multiplied by c or 1 — c, respectively. Parameters: c = 0.4, J = 0.5eV, S = §, 
W = leV , n — 0.05 

no differences in the phase diagram of the MCDA with S = \ compared to the MF 
picture but for S — § a large paramagnetic phase appears (cf. figure M again) . 
The picture will change for small J and small concentrations c. In this parameter regime 
the non-magnetic band becomes broad and the exchange splitting is so small that the 
magnetic sub-bands are lying partly or fully within the non-magnetic band. This can be 
easily seen in figure |8] with the help of the projected GFs in ([5]). In this case the effective 
filling of the magnetic sub-bands is not n eff = n/c as in the large J limit but less since 
also states of the non-magnetic band are filled. Therefore the sequence of phases with 
increasing n becomes more and more like the undiluted case at lower J. This is consistent 
because in the limit of J — > there is no difference between a diluted and a concentrated 
lattice, because both sub-systems have a vanishing potential (r] M (E) — > = ^ NM ) and 
can not be distinguished from each other. For very low couplings the same sequence 
of phases reproduces (FM: < n < 0.1, A: 0.1 < n < 0.35,C: 0.35 < n < 0.75, G: 
0.75 < n < 1) for all concentrations. Since the MF and MCDA results are the same for 
low coupling due to ( |38|) . this holds for both methods, too. 

It is worth to make some comments about the connection of a system's T = phase 
diagram and its finite temperature behavior. At large couplings J the energy differences 
between the magnetic phases are usually large (except at the borders between two 
phases). This implies that the respective phase is stable at moderate finite temperatures, 
too. On the other hand energy differences decrease when the coupling gets smaller. Thus 
even when an ordered phases appears in the phase diagram its critical temperature (Tc, 
T/v) gets very small, so that the paramagnetic phase arises at very low temperatures. 
In the extreme limit J = all phases have the same energy and are indistinguishable. 
Same is true for the limit of extreme dilution c — > 0. Even though the phase diagram 
would consist almost only of the ferromagnetic phase for n < c, what maybe contradicts 
the intuition, the Curie temperature would be very low for all couplings J, possibly even 
T c = 0+ 
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Figure 9. Phase diagrams with dilution for two different concentrations, left column: 
c = 0.4 right column: c — 0.8. The white lines show the band filling n = c. The two 
phase diagrams in the top line are calculated with the mean-field approximation. The 
others with the MCDA for spin S = | and S = §, respectively. Differences between 
MF and MCDA results appear at large couplings J and become smaller for large spins. 
Parameters: Free electron bandwidth: W = 1.0 eV. 

4. Summary and Outlook 

We have compared the internal energy of several ordered magnetic phases and the 
paramagnetic phase in a diluted ferromagnetic Kondo-lattice model. This allows us 
to make conclusions about the existence of the various magnetic phases at T = 
as a function of the coupling J, the band filling n and the concentration of the 
magnetic atoms c. To get the internal energy we applied two approximations to treat 
the ferromagnetic Kondo-lattice model. The expressions for the Green's function's 
self-energy of the (concentrated) FKLM defined effective potentials and therewith a 
dynamical alloy analogy which we included into a coherent potential approximation. 
We used this to handle the problem of the diluted disordered system. As a first 
approximation for the FKLM we used a relatively simple mean-field decoupling. But 
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already this simple ansatz gave insight in the main mechanisms of stabilizing a definite 
phase at a given parameter set. A more accurate moment conserving decoupling 
approach leads to the statement that the mean-field results are appropriate at low 
couplings or large spin quantum numbers S. For low spins the system has a more 
pronounced quantum character and we have distinct deviations from the mean-field 
picture. Both methods made clear that for very low couplings J and arbitrary 
concentrations c all magnetic phases appear at the same band fillings as in an undiluted 
system. 

It is surprising that the simple mean-field ansatz leads, under certain conditions, to 
relatively good results at T = 0. This should change when we go to finite temperatures. 
The paramagnetic phase is treated completely incorrectly in the MF approximation, 
due to the importance of fluctuations at (S z ) ^ S. Thus it is not sufficient to use 
it for, e.g., calculating Curie temperatures. To do this one has to use more accurate 
many-body methods. Expanding this theory for finite temperatures is very important to 
make predictions about possible candidates for dilute magnetic semiconductors at room 
temperature. To do this, the introduced pure KLM is surely not enough to describe 
DMS properly. For example, the magnetic and non-magnetic sites consist of different 
atoms and therefore atomic levels. This results in additional changes of the according 
density of states and will be done in forthcoming work. 

Appendix A. Coefficients and expectation values of the MCDA 

The moments of the various Green's functions X(- a = ((A ia ; c^.)) , where A ia is an 
arbitrary operator combination, are defined by 



(n-p)-fold (p)-fold 

Explicit calculation of the moments yields the coefficients of the MCDA as 
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where we suppressed the index a at every site index i for readability. We have 
abbreviated 





= (S?n la ) 
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The prefactors appearing in the equations of motion are 
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The spin's expectation values are trivial for ferromagnetic saturation and for 
paramagnetism we get them from the solution by Callen |30j . This yields 

saturation : 

((s*r) = s n 

paramagnetism : 

(St) = ((St) 3 ) = 

((s?) 2 ) =\s(s + i) 

(SrSf) =^(5 + 1) 

(SrS?S!) = -z a ((S*) 2 ) . 

The expectation values of the magnetic subsystem can be calculated via the spectral 
theorem with the according Fourier-transformed Green's functions 

(rw> = ~E / dEf.(E,fi)lmGZa(E) 

n k ^ 

(Sf niaa ) = (Sf)(n i<7 ) - f dEf4E^)lmT^ aa (E) 

71 kq J 

(Sr = -]y^E / dEf4E^)ImF™ aa (E) . 
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Note that the chemical potential /i within the Fermi functions f_(E,[i) is derived from 
the whole diluted system via the condition 



This means that only the electron number of the total system is conserved, while the 
number of electrons in the magnetic sub-system can change. 
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